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■ Context. BL Lacertae is the prototype of the blazar subclass named after it. Yet, it has occasionally shown a peculiar behaviour that has questioned 
a simple interpretation of its broad-band emission in terms of synchrotron plus synchrotron self-Compton (SSC) radiation. 
I I Aims. In the 2007-2008 observing season we carried out a new multiwavelength campaign of the Whole Earth Blazar Telescope (WEBT) on BL 
Lacertae, involving three pointings by the XMM-Newton satellite in July and December 2007, and January 2008, to study its emission properties, 
T-H ' particularly in the optical-X-ray energy range. 

^ Methods. The source was monitored in the optical-to-radio bands by 37 telescopes. The brightness level was relatively low. Some episodes of very 
T-H fast variability were detected in the optical bands. Flux changes had larger amplitude at the higher radio frequencies than at longer wavelengths. 
Results. The X-ray spectra acquired by the EPIC instrument onboard XMM-Newton are well fitted by a power law with photon index F ~ 2 
and photoelectric absorption exceeding the Galactic value. However, when taking into account the presence of a molecular cloud on the line of 
T-H sight, the EPIC data are best fitted by a double power law, implying a concave X-ray spectrum. The spectral energy distributions (SEDs) built 
with simultaneous radio-to-X-ray data at the epochs of the XMM-Newton observations suggest that the peak of the synchrotron emission lies in 
the near-IR band, and show a prominent UV excess, besides a slight soft-X-ray excess. A comparison with the SEDs corresponding to previous 
observations with X-ray satellites shows that the X-ray spectrum is very variable, since it can change from extremely steep to extremely hard, 
and can be more or less curved in intermediate states. We ascribe the UV excess to thermal emission from the accretion disc, and the other 
broad-band spectral features to the presence of two synchrotron components, with their related SSC emission. We fit the thermal emission with a 
black body law and the non-thermal components by means of a helical jet model. The fit indicates a disc temperature S 20000 K and a luminosity 
'*erg s"' . 
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1 . Introduction plasma jet closely aligned with the line of sight. The jet is gen- 
erated by a supermassive black hole fed by infall of matter from 

BL Lacertae is the prototype of one of the two blazar subclasses, accretion disc. The broad-band spectral energy distribution 

the BL Lac objects, the other subclass being represented by (ggj^^ ^ ^j^^^^ ^^^^^^ ^ ^^^^^^ ^ 

the flat-spectrum radio quasars (FSRQs). Common features of representation, shows two wide bumps. The low-energy bump, 

blazars are: 1) extreme flux variability at all wavelengths, from ^j^j^j^ ^^^^^^^ ^^^^ jj^^ ^^^-^ ^j^^ optical-UV (for some BL 

radio to r-ray frequencies, on a wide variety of time scales, ^acs up to X-ray) frequencies, is ascribed to synchrotron radi- 

from long-term (months, years) oscillations to inti-a-day vai'i- ^^-^^ relativistic electrons in the jet. The high-energy bump, 

ability (IDV); 11) high radio and optical polarization; 111) bright- ^^^^^^ng the X-ray to y-ray energies, is likely due to inverse- 

ness temperatures exceeding the Compton limit; iv) superlumi- Compton scattering of seed photons off the relativistic electrons, 

nal motion of the radio components. The commonly accepted According to the synchrotron self Compton (SSC) model, the 

paradigm foresees that their non-thermal emission comes from a ^^^^ j^^^^^^ ^1^^ synchrotron photons themselves. In con- 
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trast, the external Compton (EC) scenario foresees that seed pho- 
tons may enter the jet e ither directly from the accretion disc 



for questions regarding their availability, please contact the WEBT gio" (e-g- Sikora et al., 1994) or hot corona surrounding the disc 
President Massimo Villata (villataOoato . inaf . it). (e-g- iGhisellini & Tavecchio^ ,2009,). SSC models usually fairly 
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explain the SEDs of the low-luminosity blazars, i.e. the BL Lac 
objects, while EC models are needed to fit the SEDs of the 
FSRQs. However, recent multiwavelength studies on a num- 
ber of blazars, which included observations by the y-ray satel- 
lite AGILE, have shown that multiple SSC and/or EC compo- 
nents ar e necessary to explain the observed high-energy fluxes 
(see e.g. Chen et al., 2008;'Pucella et al.", 2008;'Vercellone et al.', 
[2009; D' Ammando et al., 2009; Donnarumma et al., 2009a). 

In addition to these two non-thermal jet components, the 
SEDs of quasar-type blazars sometimes show a "blue bump" in 
between, which is thought to be the signature of the thermal radi- 
ation emitted from the accretion disc. Indeed, the spectra of these 
objects usually display prominent broad emission lines, which 
are most likely produced by photoionization of the broad line 
region due to the disc radiation. 

On the contrary, BL Lacs are by definition almost fea- 
tureless obiects_Jec|uivalent width less than 5 A in their rest 



frame. IStickel et al ] 1199 Ih . It was hence a surprise when 



IVermeulen et al.l(ll995h . and soon after lCorbett et al.l (ll996h . dis- 



covered a broad Ho- (and H/3) emission line in the spectrum of 
BL Lacertae, whose luminosity (~ 10'*' ergs"') and full-width 
half-maximum (~ 4000 km s"') are comparable to those of type 
I Seyfert galaxies such as NG C 4151. Subsequent spectroscopic 
monitoring of this source bv ICorbett et"an (l2000h showed that 
the Ha equivalent width is approximately inversely proportional 
to the optical continuum flux. This suggested that the broad line 
region is photoionized by a radiation source that is not the same 
producing the optical continuum. The photoionising radiation 
would most likely come from the accretion disc. 

Another importan t issue comes from the results obtained by 
iRavasio et al.l (l2003h . They analyzed the X-ray data acquired 
by BeppoSAX from October 31 to November 2, 2000, during 
an extensive multiwavelength campaign. When constructing the 
source SED with contemporaneous data, it was evident that the 
steep X-ray spectrum was offset with respect to the extrapolation 
of the optical one. One possible explanation was the presence of 
an extra component in addition to the synchrotron and inverse- 
Compton ones. 

Finally, acco rding to iMadeisld et al.| (1 19991) 
iBottcher & BloomI dlOOOt) . the explanation of the 



and 
y-ray 



flux detected by the EGRET instrument on boa rd the CGRO 
satell ite during the 1997 optical outburst (see iBloom et al.L 
1 19971) requires an EC emission component in addition to the 
SSC one. 

Taken together, these results suggest that in BL Lacertae the 
interpretation of the broad-band emission may require a more 
complex scenario than that usually envisaged for the BL Lac 
objects, involving just one synchrotron plus its SSC emission. 

In the last decade, BL Lacertae has been extensively stud- 
ied by the Whole Earth Blazar Telescope (WEBT) collabo- 
ratiorOl which has carried out several multiwavelength cam- 
3aigns o n this object ( Villata ej al., 2002; Rayasio et al., 2002; 




2003; Villat a et all 12004, 
2007^ .ViUata et all 1200' 



at [Bach et all 120061; 
collecting tens of 



thousands of optical-to-radio data. These studies were focused 
on its multiwavelength flux variability, colour behaviour, corre- 
lations among flux variations in different bands, possible peri- 
odicity of the radio outbursts. The main aim of the new WEBT 
campaign organized in the 2007-2008 observing season was in- 
stead to address the problem of disentangling the possible mul- 
tiple contributions to the BL Lac flux from the radio band to y- 
rays. For this sake, the optical-to-radio monitoring by the WEBT 
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was complemented by three pointings by the XMM-Newton 
satellite. Moreover, we also obtained optical spectra with the 3.6 
m Telescopio Nazionale Galileo (TNG) to investigate the prop- 
erties of the Ha broad emission line, and possibly infer informa- 
tion on the accretion disc. The results of the spectroscopic study 
will be reported elsewhere. 

This paper is organised as follows. In Sect. 2 we present the 
WEBT optical-to-radio light curves. The analysis of the XMM- 
Newton data is reported in Sect. 3. In Sect. 4 we show the SEDs 
corresponding to the XMM-Newton epochs, and compare them 
with those related to previous observations by X-ray satellites. 
The interpretation of the XMM-Newton SEDs is discussed in 
Sect. 5. Finally, Sect. 6 contains a summary and discussion of 
the main results. 



2. Multifrequency observations by the WEBT 

The new WEBT campaign on BL Lacertae took place in the 
2007-2008 observing season. The participating observatories 
are listed in Table [T] Optical and near-IR data were collected 
as instrumental magnitudes of the source and reference stars 
in the same field to apply th e same calibration dBertaud et al.L 
1969; Fi orucci & Tostilll996 1. The light curves obtained by as- 
sembling all datasets were carefully inspected to correct for sys- 
tematic off'sets and to reduce data scatter by binning noisy data 
taken by the same observer within a few minutes. The results 
are shown in Fig.[T] where the vertical lines indicate the epochs 
of the three XMM-Newton pointings. We can see a noticeable 
flux variability, which progressively increases its amplitude go- 
ing from the / to the B band. By considering only the period 
of common monitoring (before JD - 2454500), the maximum 
variability amplitude (maximum - minimum) is 1 .42, 1 .46, 1 .50, 
and 1.55 mag in the /, R, V, and B bands, respectively. The near- 
IR time coverage is inferior to the optical one, but the near-IR 
data are important to add information to the SED (see Sect. 5). 

The majority of the variability episodes have a time scale of 
a few days, but we can also recognise a long-term increasing 
trend starting from JD ~ 2454300, as well as a few very fast 
events. One of these involved a brightening of about 0.9 mag in 
24 hours, from R ~ 14.6 on JD = 2454301.5 to R ~ 13.7 the 
night after, when observations at the Valle d'Aosta Observatory 
showed a source brightening of ~ 0.3 mag in less than 3 hours. 
This behaviour was confirmed by observations in the V and 
/ bands, ruling out that this rapid flux increase was an arti- 
fact. Similar fast variations are not uncommon in BL Lacertae. 
Wh en analys i ng the 13248 j^-ban d data acquired by the WEBT 
(Villata et all [200ll2004blaL [2009. and this paper) during 1500 
nights over more than 15 years, though with inhomogeneous 
sampling, we can distinguish between two kinds of rapid flux 
variability. 

- Fast and noticeable intraday variations: we consider varia- 
tions > 0.25 mag with rate > 0.1 mag/hour. These were 
found in 25 nights out of 677 nights where the observing 
time coverage is > 2.5 hours, with the maximum amplitude 
episode involving a change AR - 0.52 in about 3.9 hours. 

- Large interday variations: we consider variations > 0.75 mag 
in < 36 hours. These were observed 6 times, and the above 
mentioned episode (about 0.9 mag in 24 hours) is the most 
extreme one. 

Radio data were collected as already calibrated flux den- 
sities. The radio light curves are shown in Fig. [2l where the 
first panel displays the /?-band light curve for comparison. We 
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Fig. 1. Optical UBVRI and near-IR JHK light 
curves of BL Lacertae in the 2007-2008 ob- 
serving season. Vertical lines indicate the three 
XMM-Newton pointings of July 10-11 and 
December 5, 2007, and January 8, 2008. 



also included data from the VLA/VLBA Polarization Calibration 
Database (PCDfl As expected, the radio flux variations, which 
are more evident at the shortest wavelengths, are smoother than 
the optical variation^, and the radio time scales are longer. 
Moreover, the long-term increasing trend characterising the op- 
ti cal light curve is not recognisable in the radio band. According 
to lVillata et al.l (|2009|) . the optical outbursts of BL Lacertae are 
usually followed by high-frequency radio events, with time de- 
lays of at least 100 days, which can grow to 200 or even 300 
days, depending on the relative orientation of the coiTespond- 
ing emitting regions in the jet. Hence, the high optical level ob- 
served at the beginning of our observing period, in May-June 
2007, might be related to the bright radio state that is visible at 
the higher radio frequencies around JD = 245440(fl 



^ |http : //www . via . nrao . edu/astro/calib/polar/ 1 

^ This is even more evident when comparing radio flux densities to 
optical flux densities instead of magnitudes. 

* Indeed, the PCD shows no further radio event until May 2009, the 
43 GHz flux density remaining below 3 Jy. 



3. Observations by XMM-Newton 

The X-ray Multi-Mirror Mission (XMM) - Newton satellite ob- 
served the source on July 10-1 1 and December 5, 2007, and then 
on January 8, 2008 (PI: C. M. Raiteri). 



3.1. EPIC data 

The European Photon Imaging Camera (EPIC) onb oard XMM- 
Newton includ es three detec t ors: M OSl, MOS2 ( T urner et all 
2001), and pn (IStriider et al.L |2001|) . Since a bright state of the 
source could not be excluded, we chose a medium filter to 
avoid possible contamination by lower-energy photons; more- 
over, we selected a small-window configuration to minimize 
possible photon pile-up. 

Data were reduced with the Science Analysis System (SAS) 
software, v ersion 8.0.1, f ollowin g the same standard procedure 
adopted in 'R aiteri et al.l (12007b'). A temporal filtering was ap- 
plied to discard high-background periods. Source spectra were 
extracted from circular regions of ~ 40 arcsec radius. The MOS 
background was estimated on external CCDs, while for the pn 
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Fig. 2. i?-band liglit curve of BL Lacertae in 
2007-2008 (top panel) compared to tlie be- 
haviour of tire radio flux densities (Jy) at dif- 
ferent frequencies. Vertical lines indicate the 
three XMM-Newton pointings of July 10-11 
and December 5, 2007, and January 8, 2008. 



we selected the largest source-free circle that could be arranged 
on the same CCD. 

The source spectra were binned with the grppha task of the 
FTOOL package, to have a minimum of 25 counts in each bin. 
The binned spectra were then analysed in the 0.35-12 keV en- 
ergy range with the Xspec task of the XANADU package. We 
fitted the MOSl, MOS2, and pn spectra of each epoch together 
to increase the statistics. 

We first considered a single power law model with free ab- 
sorption, where the absorption is modelled according to the 
[Wilms et al. (2000) prescriptions. The results of this spectral 
fitting for the three epochs are displayed in the top panels of 
Figs. [3}|5] The corresponding model parameters are reported in 
Table|2] where Col. 2 gives the hydrogen column density. Col. 3 
the photon spectral index F, Col. 4 the unabsorbed flux density at 
1 keV, Col. 5 the 2-10 keV observed flux, and Col. 6 the value of 
X^/v (being v the number of degrees of freedom ). The ;\f^/v val- 
ues in Table |2] indicate that the model is acceptable. This is also 
confirmed by the deviations of the observed data from the folded 
model, which are plotted in the bottom panels of Figs.[3]-l5] The 
best-fit value of the hydrogen column density varies between 2.8 
and 3.1x10^' cm"^, within the range of values found in previous 
analyses (see Sect. 5). 

Actually, the Galactic atomic hydrogen column density 
toward BL Lacertae is A^h - 1.71 x 1 0^^ cm~^ (from 
the Leiden/ Argentine/Bonn (LAB) Survey, see iKalberla et al.L 



'2005*). However, observations of local interste llar CO to- 
ward BL Lac have revealed a molecular cloud ( iBania et all 
, 1991l;lLucas & Liszill993l;lLiszt & Lucaslll998l) . According to 
Liszt &Lucasl (1 199 81). its '^CO column density is (8.48 + 0.78) x 



lO'^^cm""^. Assuming that the molecular hydrog en column den- 
sity N h^ is typically 10^ times the '^CO one dLucas & Liszd 
1 19931) . we derive an hydrogen column density of ~ 1.7 x 
10^' cm"^ due to the molecular cloud. This value depends on 
the uncertain ratio between CO and h£], but taking it at face 
value, the total hydrogen column density toward BL Lac be- 
comes A^H = 3.4 X 10^' cm"^. This value is not very far from the 
A^H values we found when fitting the EPIC spectra with a power 
law model with free absorption, but as we will see in Sect. 5 this 
modest difference can make a difference in the interpretation of 
the source X-ray spectrum. 

We thus fixed A^h - 3.4 x 10^' cm"^ and re-fitted a single 
power law model to the EPIC spectra. As expected, the goodness 
of the new fits is inferior to the previous case, and a slight excess 
of counts in the soft X-ray domain appears. This could be the 
signature of a curvature in the source spectra. Indeed, when we 
adopt a double power law model, the fit improves significantly, 
as is shown in Table [3] and this suggests that the spectrum is 
concave. 



From the recent paper bv lLisztl ( l2007l) we can infer that A'hj /Nij 
is most likely in the range ~ 1-2 x 10*. 
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Table 1. List of optical, near-IR, and radio observatories con- 
tributing data to this work. 



M0S1+M0S2+pn 



XMM-Newton: July 1 0-1 1 , 2007 



Observatory 


Tel. size 


Bands 




Optical 




Abastumani, Georgia 


70 cm 


R 


ARIES, India 


104 cm 


BVRI 


Armenzano, Italy 


35 cm 


BRI 


Armenzano, Italy 


40 cm 


BVRI 


Belogradchik, Bulgaria 


60 cm 


VRI 


BOOTES-2, Spain 


30 cm 


R 


Bordeaux, France 


20 cm 


V 


Calai- Alto, Spain" 


220 cm 


R 


Crimean, Ukraine 


70 cm 


BVRI 


Kitt Peak (MDM), USA 


130 cm 


UBVRI 


L'AmpoUa, Spain 


36 cm 


R 


Lulin (SLT), Taiwan 


40 cm 


R 


Michael Adrian, Germany 


120 cm 


R 


Mt. Lemmon, USA 


100 cm 


BVRI 


New Mexico Skies, USA 


30 cm 


VRI 


Roque (KVA), Spain 


35 cm 


R 


Roque (NOT), Spain 


256 cm 


UBVRI 


Rozhen, Bulgaria 


50/70 cm 


BVR 


Rozhen, Bulgaria 


200 cm 


BVRI 


Sabadell, Spain 


50 cm 


R 


San Pedro Martir, Mexico 


84 cm 


R 


Sobaeksan, South Korea 


61 cm 


VRI 


Sommers-Bausch, USA 


61 cm 


VRI 


St. Petersburg, Russia 


40 cm 


BVRI 


Talmassons, Italy 


35 cm 


BVR 


Teide (BRT), Spain 


35 cm 


BVR 


Torino, Italy 


105 cm 


BVRI 


Tuorla, Finland 


103 cm 


R 


Univ. of Victoria, Canada 


50 cm 


R 


Valle d'Aosta, Italy 


81 cm 


BVRI 




Near-infrared 




Campo Imperatore, Italy 


110 cm 


JHK 


Roque (NOT), Spain 


256 cm 


HK 




Radio 




Crimean (RT-22), Ukraine 


22 m 


37 GHz 


Mauna Kea (SMA), USA 


8 x6m'' 


230, 345 GHz 


Medicina, Italy 


32 m 


5, 8, 22 GHz 


Metsahovi, Finland 


14 m 


37 GHz 


Noto, Italy 


32 m 


43 GHz 


UMRAO, USA 


26 m 


4.8, 8.0, 14.5 GHz 



■ Calar Alto data were acquired as part of the MAPCAT (Monitoring 
AGN with Polarimetry at the Calar Alto Telescopes) project. 
* Radio interferometer including 8 dishes of 6 m size. 



The Iv values in Table [3] are a bit smaller than those in 
Table 111 To better compare the two model fits, we calculated the 
F-test probability, which is 1.25x10-^ for July 10-1 1, 2.70x10-^ 
for December 5, and 3.11 x 10""' for January 8, 2008. These 
results suggest that the double power law model with fixed ab- 
sorption may be more appropriate to describe the EPIC spectra 
than the single power law model with free absorption. 



3.2. OM data 

XMM-Newton also ca rries an optical-UV 30 cm telescope, the 
Optical Monitor (OM, Mason et an.l200lh . The BL Lac observa- 
tions were performed with all its filters: V, B, U, UVWl, UVM2, 
and UVW2, with long exposures (see Table |4]i. The OM data 
were reduced with the SAS software, version 8.0.1. The tasks 
omsource and omphotom were used to perform aperture pho- 
tometry on the images produced by omi chain. The resulting 
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Fig. 3. EPIC spectrum of BL Lacertae on July 10-11, 2007; 
black squares, red triangles, and green diamonds represent 
MOSl, M0S2, and pn data, respectively. The bottom panel 
shows the deviations of the observed data from the folded model 
(a power law with free absorption) in unit of standard deviations. 
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XMM-Newton: December 5, 2007 
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Fig. 4. EPIC spectrum of BL Lacertae on December 5, 2007; 
black squares, red triangles, and green diamonds represent 
MOSl, M0S2, and pn data, respectively. The bottom panel 
shows the deviations of the observed data from the folded model 
(a power law with free absorption) in unit of standard deviations. 

magnitudes are reported in Table |4] The uncertainties take into 
account the measure, systematic and calibration eiTors. In the 
optical filters, where a comparison with ground-based measure- 
ments is possible, the OM magnitudes of the reference stars (B 
C H K) are within 0. 1 mag with respect to the values we adopted 
for the calibration of the ground data, but they are stable (within 
2-3 hundredths of mag) in the three XMM-Newton epochs. 

4. Spectral energy distributions 

The upper panel of Fig. |6] displays the broad-band SEDs cor- 
responding to the XMM-Newton observations analysed in the 
previous section. We show both the single power law with free 
absorption and the double power law with atomic plus molecu- 
lar Galactic absorption fits to the X-ray spectra. Optical and UV 
magnitudes were coiTected for Galactic extinction by adopting 
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Table 2. Results of fitting the EPIC data with a single power law with free absorption. 



Date 


Na 


r 


^IkeV 










[10" cm~2] 






[erg cm~^ s 






(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


July 10-11, 2007 


3.05 ± 0.06 


2.01 + 0.02 


2.58 ± 0.04 


9.64 X 10- 


12 


0.942 (1347) 


December 5, 2007 


2.92 + 0.06 


1.99 + 0.02 


1.96 + 0.04 


7.65 X 10" 


12 


1.026 (1250) 


January 8, 2008 


2.86 ± 0.06 


1.91 ±0.01 


1.95 + 0.03 


8.51 X 10- 


12 


0.905 (1381) 



Table 3. Results of fitting the EPIC data with a double power law with fixed total atomic and molecular column density. The two 
power laws have photon indices and and unabsorbed flux densities at 1 keV F^^^^^ and F^-^.^^, respectively 



Date 


A'h 


F' 


F^ 


IkcV 


IkcV 


-P'2-iokev 






[10" cm-2] 










[ergcm^^ s"'] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 






July 10-11, 2007 


3.40 


2.48 


1.72 


1.63 


1.14 


9.76 X 10"'^ 


0.939 (1346) 


December 5, 2007 


3.40 


2.58 


1.67 


1.26 


0.90 


7.77 X 10-'2 


1.013 (1249) 


January 8, 2008 


3.40 


2.48 


1.51 


1.45 


0.73 


8.69 X 10-'2 


0.880(1380) 



Table 4. Optical-UV magnitudes of BL Lacertae derived from the data analysis of the OM frames. 



Date 


V 


B 


U 


UVIVI 


UVM2 


UV1V2 






Exposure times (s) 








July 10-11, 2007 


1498 


1499 


1499 


2601 


3098 


6257" 


December 5, 2007 


2100 


2099 


2099 


3500 


3780 


4000 


January 8, 2008 


1700 


1700 


1698 


2799 


3300 


7359* 








Magnitudes 








July 10-11, 2007 


15.33 + 0.10 


16.27 + 0.10 


15.87 + 0.10 


15.99 + 0.10 


16.90 + 0.10 


17.13 + 0.14 


December 5, 2007 


14.96 + 0.10 


15.91 +0.10 


15.52 + 0.10 


15.62 + 0.10 


16.45 + 0.10 


16.98 + 0.11 


January 8, 2008 


14.97 + 0.10 


15.90 + 0.10 


15.53 + 0.10 


15.65 + 0.10 


16.47 + 0.10 


16.88 + 0.14 



" Two exposures of 3579 and 2678 s. 
* Two exposures of 3779 and 3580 s. 



XMM-Newton: January 8, 2008 
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Fig. 5. EPIC spectrum of BL Lacertae on January 8, 2008; black 
squares, red triangles, and green diamonds represent MOSl, 
MOS2, and pn data, respectively. The bottom panel shows the 
deviations of the observed data from the folded model (a power 
law with free absorption) in unit of standard deviations. 



Ab - 1 .4-2 from lSchlegel et akl (Il998l) and calculating the v alues 
at the other wavelengths according to lCardelli et alj(ll989l) . De- 
reddened ground-based optical magnitudes were then converted 
into fl uxes using the zero-mag fluxes given by BesselLet al. 
( 119981) : as for the optical and UV magnitudes from the OM, this 
conversion was performed using Vega as calibrator 



The optical flux densities were further coiTected for the con- 
tribution of the host galaxy. Assuming an /?-band magnitude of 
15.55 after [Scarpa e t al.' (2000) and the average colour indices 
for elliptical galaxies by Mannucci et al. (2001), the host galaxy 
flux densities are: 10.62, 13.97, 11.83, 5.90, 4.23, 2.89, 1.30, and 
0.36 mJy in K, H, J, I, R, V, B, and U bands, respectively. Using 
a De Vaucouleurs' profile, we estimated that the contribution to 
the observed fluxes is ~ 60% of the whole galaxy flux (see also 
Villata et al., 2002). In the UV, the results of spectral evolution 
modelling of stellar populations bv lBruzual & Chariot! (l2003b al- 
lowed us to estimate that the contribution from the host galaxy 
may be neglected. Indeed, if we consider ages between 4 and 
13 Gyr, the galaxy flux density at ~ 2000 A is about 50 to 100 
times lower than in the R band, i.e. the contribution of the galaxy 
would affect the measured flux in the UVW2 band by ~ 1-2%. 

Optical data are strictly simultaneous to the XMM-Newton 
observations. Most radio data are simultaneous too, but in some 
cases we considered data taken a few days earlier or later. We 
note the excellent agreement between the V, B, and U data taken 
with ground-based telescopes and the corresponding data ac- 
quired by the OM. The main features of the SEDs in the figure 
are: 

- the peak of the synchrotron emission lies in the infrared (see 
below); 

- since a higher radio brightness coiTesponds to a lower optical 
state and viceversa, the synchrotron peak likely shifts toward 
higher frequencies as the optical flux increases; 

- there is a strong UV excess, since the UV points do not lie 
on the extrapolation of the optical trend; 
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Fig. 6. Top panel: broad-band SEDs of BL 
Lacertae during the three XMM-Newton 
observations of July and December 2007, and 
January 2008. Both the single power law with 
free absorption and the double power law with 
atomic plus molecular Galactic absorption 
fits to the X-ray spectra are shown. Bottom 
panel: SEDs corresponding to epochs of pre- 
vious X-ray satellites observations: ASCA i n 
November 1995 (blue, Sambruna et al.','l999!); 
RXT E in July 1997 (grey, Madejski et aU 
1999); ASCA -l- RXTE in July 1997 (red, 
Tanihata et al.L l2000h : B eppoSAX in 
Nove mber 1997 (light green, jPadovani et all 

1200 Ih : BeppoSAX in June (cyan ) and 

December (pink) 1999 jRavasio et all l2002h : 
BeppoSAX in July (orange) and November 
(green) 2000 jRavasio et all l2003h . All radio 
and optical data are from the WEBT archive; 
the vertical bars on the red optical points 
indicate the range of variation of the optical 
fluxes during the ASCA-l-RXTE observations 
in July 1997. The r ange of flux densitie s 
detected by EGRET jlJartman et all 1 19991) . 
as well as the TeV spect rum observed by 
MAGIC jAlbertet all 120071) are also reported. 
The solid lines represent the sum of different 
emission contributions: a low-energy syn- 
chrotron-l-SSC component (dotted lines), a 
high-energy synchrotron+SSC component 
(dashed lines), and a thermal component from 
an accretion disc (dotted-dashed line) with a 
temperature of ~ 20000 K and a luminosity of 
6 X lO"*"* erg s"' (see the text for details). 



- the X-ray spectrum has either a null slope, or, more likely, it 
is concave, producing a mild soft-X-ray excess, and suggest- 
ing that in this energy region two emission components are 
intersecting each other; 

- the optical steepness makes an interpretation of the soft X- 
ray excess in terms of the tail of the synchrotron component 
unlikely. 

Details of the near-IR-to-UV SED of BL Lacertae are shown 
in Fig.|2l Here we can better appreciate the features of the three 
SEDs obtained at the epochs of the XMM-Newton observations. 
Moreover, we also report the total uncertainty affecting the UV 
points, when considering both the error on the data and the sam- 
ple variance about th e Galac tic mean extinction curve, following 
iFitzpatrick & Mas"sal (l2007l) . As stressed by these authors, the 
sample variance must be taken into account if we want to esti- 
mate a realistic error on the de-reddened SEDs. In our case, the 
lower limit to the UV fluxes that we obtain when considering the 
total uncertainty indicates that the UV excess might be smaller 
than shown by the points, but it exists, since the UV points can- 
not be shifted down enough to lie on the extrapolation of the 
optical, (quasi-power law) synchrotron trend. 

Figure |7] also shows two other SEDs built with simultane- 
ous near-IR and optical data. The near-IR data have been de- 



reddened and cleaned from the host galaxy contribution simi- 
larly to the optical data. In the optical frequency range, these 
two SEDs confirm the trend shown by the SEDs obtained at the 
three XMM-Newton epochs, but they add important information 
in the near-IR. Indeed, they suggest that the peak of the syn- 
chrotron component lies in this energy range, or close by. 

4.1. Comparison with previous observations 

The bottom panel of Fig.|6]shows SEDs coiTesponding to epochs 
when various X-ray satellites observed BL Lacertae. The X-ray 
spectral fits were taken from the literature, while we searched the 
massive WEBT archive on this source, containin g all radio-to- 
optical data f r om the previous WEBT campaigns jVillata et all 
120021 12004blli 120091) as well as data from the literature, for con- 
temporaneous low-energy data. 

When considering the behaviour of the X-ray spectrum, it is 
not easy to compare our results to those found by other authors 
when analysing different X-ray data. The reason is that the shape 
of the X-ray spectrum strongly depends on the choice of the ab- 
sorption, i.e. on the A^h value adopted to perform the spectral 
fits. Models where the hydrogen column density is left to vary 
freely yield a variety of A^h values (mostly in the range 1 .4-3.5 x 
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Fig. 7. Details of the near-IR-to-UV SED of BL Lacertae. In the 
three optical-UV SEDs obtained at the epochs of the XMM- 
Newton observations of July 10-11, 2007, December 5, 2007, 
and January 8, 2008, empty symbols represent ground-based ob- 
servations, while filled symbols refer to data taken by the OM. 
The dotted error bars on the UV points indicate the total un- 
certainty, resulting from considering both the data error (solid 
error bars) and the sample variance about the mean extinction 
curve, according to lFitzpatrick & Massal(l2007l) . To clearly show 
the uncertainties holding on the UV data, the UV points corre- 
sponding to the December 2007 and January 2008 epochs have 
been slightly shifted in frequency around the effective value. For 
comparison, we show other two SEDs obtained with data ac- 
quired by ground-based telescopes on June 24 and October 29- 
30, 2007, when near-IR data contemporaneous to the optical data 
were available. The near-IR portion of the SEDs shows that the 
peak of the synchrotron emission probably lies in this energy 
range, or close by. 



10^' cm^^). These are usually lower than the estimated Galactic 
total absorption, due to both atomic and molecular hydrogen. 
The point is that the amount of absorption due to molecular hy- 
drogen is not directly measurable, whi ch introduces a further un- 
certainty (see Sect. 4 and discussion in lMadejski et al.Lll999h . In 
previous analyses of X-ray data, different authors have followed 
different prescriptions. An A^h value close to the one we assumed 
was also investigated bv lSambruna et al] d 19991) when analysing 
ASCA observations in November 1995, and bv iPadovani et al.l 
(1200 lb for observations performed by BeppoSAX in November 
1997 (and by ROSA T in 1992). In b oth epochs the X-ray flux 
was relatively low and lSambruna et aL (1999) found that the best 
fit model was a broken power law, leading to a concave X-ray 
spectrum. 

Observations by RXTE and A SCA during the big out- 
burst of July 1997 w ere analysed by Madeiski et alj d 19991) and 
iTanihata et all (|2000|) . The latter authors combined the data of 
both sateUites and distinguished between a low and a high state. 
They used a high A^h - 4.6 x 10^' cm"^, and found that the high 
state was best fitted by a double power law model with an ex- 
tremely stee p spectral index bel ow ~ 1 keV. In contrast to this 
high value, iRavasio et al.1 (l2002h and iRavasio et al.l (l2003l) pre- 
ferred an A^H - 2.5 x 10^' cm"^ to analyse BeppoSAX obser- 
vations in 1999-2000. Nevertheless, the June 1999 data were 
best fitted by a double power law model, implying a notice- 
able upward curvature in the soft X-ray domain. In contrast, 
in November 2000 the whole X-ray spectrum was steep, and 



its extrapolation to optical frequencies was not intersecting the 
simultaneous optical spectrum. This offset was questioned by 
Bottcher et al. (2003), who warned about the possibility that 
it was produced by flux averaging in a period of noticeable 
variability, as shown by the intensive monitoring performed by 
the WEBT observe rs during the 2000-2001 observing season 
(IVillata et all 12002 *). However, to remove the offset we should 
have missed one or more big flares, such as to increase the mean 
optical level by more than 1 mag, which seems unlikely. Hence, 
if in previous works a soft spectral component at low X-ray en- 
ergies had always been interpreted as the tail of the synchrotron 
emission, this was no longer a plausible explanation for the 
November 2000 SED, as it is not a plausible explanation for the 
XMM-Newton SE Ds presented above. A number of hypotheses 
were suggested by IRavasio et al.l ( 1200 3^ to justify the optical- 
X-ray offset: a sudden increase of the dust-to-gas ratio toward 
BL Lacertae; the detection of a bulk Compton emission; the in- 
terplay of two different synchrotron components; Klein-Nishina 
effect on the synchrotron spectrum. 

5. The helical jet model 

The availability of broad-band SEDs built with simultaneous 
data, including UV information, at the epochs of the XMM- 
Newton observations of 2007-2008 is a formidable tool to in- 
vestigate the nature of BL Lac multiwavelength emission. The 
picture that we described in the previous section suggests that the 
SED behaviour cannot be explained in terms of one synchrotron 
plus its SSC emission components. 

Thermal emission from an accretion disc could account 
for the UV excess. Indeed, the optical spectrum of this 
source occasionally shows broad emission lines that are pro- 
duced in the broad line region, which is most likely pho- 
toionised b y the radiation coming from the accretion disc 
dVermeulen et al.l 119951; fcorbett et al.L 1199a 12000) . Moreover, 
a big blue bump produced in the SED by thermal emission from 
an accretion disc has a lready been found for other quasar-type 
blazars, like 3C 273 (Smith et al., '1993"; von Montign v et all 
1997; Grandi & Palumbo, 2004; Turler et al.^2006),' 3C 279 
(Plan et al., 1999), 3C 345 ( Bregmanet al.[ll986h . and 3C 454.3 
(.Raiteri et aL .2007b. .2008.) . 

However, a thermal emission component would hardly be 
able to explain the extremely variable X-ray spectrum shown in 
the bottom panel of Fig. |6] Indeed, at 0.3 keV BeppoSAX ob- 
served a flux variation of a factor ~ 50 in about 1 year We notice 
that the behaviour of BL Lacerta e in the U V-X-ray en ergy range 
is similar to that found by RaiterTetaD d2005L |2006b a) when 
analysing the SEDs of another BL Lac object, AO 0235+164, 
which also occasion ally shows broad emission lines (see e.g. 
iRaiteri et al.L 1200731 and references therein). In that case, the al- 
ternative hypothesis of a second, higher-frequency synchrotron 
component, likely coming from an inner region of the jet, was 
left open. 

We investigated thi s problem by mean s of the helical jet 
model by Villata & Rai teril ([T9991 see also iRaiteri et al.llT999l 
IRaiteri & Vill ata 2003, and' Ostorero et al.ll2004 . We tried to re- 
produce the "XMM-Newton" broad-band SEDs of December 
2007 and January 2008, which are very similar To better trace 
the SED shape, we also added synthetic near-IR data, derived 
by shifting the October 2007 SED of Fig. [T] to match the above 
"XMM-Newton" SEDs. This implied a shift of log(vFy) - 
-0.07. 

The helical jet model presented in IVillata & Raiteril d 19991) 
foresees that orbital motion in a binary black hole system, cou- 
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pled with the interaction of the plasma jet with the surrounding 
medium, twists the emitting jet in a rotating helical structure. 
We here recall the main features of the model. The axis of the 
helical-shaped jet is assumed to lie along the z-axis of a 3-D ref- 
erence frame. The pitch angle is ^ and \j/ is the angle defined by 
the helix axis with the line of sight. The non-dimensional length 
of the helical path can be expressed in terms of the z coordinate 
along the helix axis: 



Kz) 



cos ^ 



0<z< 1, 



(1) 



which corresponds to an azimuthal angle if{z) - az, where the 
angle a is a constant. The jet viewing angle varies along the he- 
lical path as 



cos 0{z) - COS ip COS ^ -H sin (/r sin ^cos(0 - az) , 



(2) 



where (p is the azimuthal difference between the line of sight and 
the initial direction of the helical path. 

The jet is inhomogeneous: it emits radiation at progressively 
increasing wavelengths by proceeding from its apex outwards. 
Each slice of the jet can radiate, in the plasma rest reference 
frame, synchrotron photons from a minimum frequency to 
a maximum one v^^^^. Both these frequency limits decrease for 
increasing distance from the jet apex following a power law: 



1 + 



>o. 



(3) 



where l, are length scales, and / = min, max. The high-energy 
emission producing the second bump in the blazars SED is as- 
sumed to be the result of inverse Compton scattering of the syn- 
chrotron photons by the same relativistic electrons emitting them 
(SSC model). Consequently, each portion of the jet emitting syn- 
chrotron radiation between v'^^^J^l) and v'^^,^^{l) will also pro- 
duce inverse Compton radiation between n)in(0 and maxW^ 
with .{l) — \'y](J)v[ fl). The electron Lorentz factor ranges 
from yniin = 1 to yniax(^)> which has a similar power law depen- 
dence as in Eq.[3] with power Cy and length scale ly. As photon 
energies increase, the classical Thomson scattering cross section 
is replaced by the Klein-Nishina one, which takes into account 
quantum effects. Its consequence is to reduce the cross section 
from its classical value, so that Compton scattering becomes less 
efficient at high energies. We approximated this effect by requir- 
ing that V^^^.^J^l) be averaged with v'^^.^^{r) = '-^ rmax(0 when 

rmax(Ov;,max(0 > 

We assume a power law dependence of the observed flux 
density on the frequency and a cubic dependence on the Doppler 
beaming factor 6: Fy{v) oc S^v^"", where oq is the power law 
index of the local synchrotron spectrum, 6 = [r( 1 -y6 cos ff)]^^,/3 
is the bulk velocity of the emitting plasma in units of the speed of 
light, r = (1 -ffy^^^ the corresponding bulk Lorentz factor, and 
6 is the viewing angle of Eq. |2l Since the viewing angle varies 
along the helical path, also the beaming factor does. Hence, the 
flux at V peaks when the part of the jet mostly contributing to it 
has minimum 0. 

The emissivity decreases along the jet: both the synchrotron 
and inverse Compton flux densities are allowed to drop when 
moving from the jet apex outwards. For a jet slice of thickness 
dl: 



di^v,s(v) (X 6\l) V' 



I 

1 + - 

h 



dl, o , 



(4) 



Table 5. Main parameters of the helical model for the fit to the 
"XMM-Newton" broad-band SED of December 2007 - January 
2008. The differences between the low- and high-energy syn- 
chrotron+SSC components are easily seen. 



Parameter 


Low 


High 




30° 


30° 




25° 


25° 


a 


110° 


110° 


<P 


-8° 


20° 


log<(0) 


14.0 


17.8 


^ iTiin,max 


2.5 


2.5 


log /rain 


-3.2 


-3.2 


log /raax 


-1.6 


-1.6 


logrmax(O) 


3.5 


4.4 


Cy 


1.25 


1.25 


lOg/y 


-1.6 


-1.6 




0.5 


0.5 


r 


10 


10 


Cs 


1 


1 


log/s 


-1 


-1 


Cc 


1 


1 


log/. 


-1 


-1 



dF,,e(v)<x <J^(Oy-"° (1 + - 



In 



x(0 



(0 



dl. 



Cc > (5) 



For both the synchrotron and inverse Compton components, 
the observed flux densities at frequency v coming from the 
whole jet are obtained by integrating over all the jet portions 
Az,(v) contributing to that observed frequency, i.e. for which 
6{z)v'^^^{z) < V < 6(z)v'^^^{z). The total observed flux density at 
frequency v is finally obtained by summing the synchrotron and 
inverse Compton contributions. Notice that the intrinsic jet emis- 
sion does not vary with time, but the observed one may change 
as the orientation changes. 

The fit to the "XMM-Newton" broad-band SED of 
December 2007 - January 2008 in Fig.|6](solid line) is obtained 
by considering two synchrotron emission components from dif- 
ferent regions of a helical jet, with their corresponding SSC, plus 
a thermal component, modelled as a black body. The main pa- 
rameters of the model are reported in table |5] 

The lower-energy synchrotron+SSC emission (dotted line) 
comes from a helical portion that is initially fairly aligned with 
the fine of sight (4> - -8°, 0(0) ~ 6.2°), while the higher-energy 
synchrotron+SSC emission (dashed line) is produced by another 
helical region that is initially less aligned with the line of sight 
(<A = 20°, 0(0) ^ 10.4°). 

A thermal component that fits the UV excess must have a 
black body temperature ^ 20000 K and a luminosity ^ 6 x 
10"''* erg s"'. This lower limit to the temperature (and conse- 
quently to the luminosity) is constrained by the break of the 
SED in the optical-UV transition, but much hotter and hence 
more lumino us discs are pos sible. For comparison, the thermal 
disc fitted bv iPian et al.l (Il999.) to the UV data of 3C 279 has a 
temperature of 20000 K and a luminosity of 2 x 10''^ erg s ' . 

The two synchrotron+SSC components can change a lot for 
variations e.g. of the angle <p, which happens if the helix rotates, 
allowing us to explain the noticeable spectral variability of the 
source also in the absence of intrinsic, energetic processes. With 
reference to the SEDs shown in the bottom panel of Fig. |6] the 
high-energy synchrotron+SSC emission was giving an excep- 
tional contribution in November 2000, while it was very faint in 
December 1999. A detailed investigation of the model parame- 
ter space to fit the SED shape of BL Lacertae at different epochs 
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goes beyond the scope of this paper. We notice however that our 
model fit produces a GeV spectrum with photon index F ^ 2 and 
can fairly reprodu ce the TeV sp ectrum observed by the MAGIC 
telescope in 2005 d Albert et al.Ll2007h . 

6. Discussion and conclusions 

The WEBT campaign on BL Lacertae in the 2007-2008 ob- 
serving season involved 37 optical-to-radio telescopes. They ob- 
served the source in a relatively faint state. Nevertheless, some 
fast variability episodes were detected in the optical bands, su- 
perposed to a long-term flux increasing trend. During the cam- 
paign, three observations by the XMM-Newton satellite added 
information on the UV and X-ray states of the source. 

The broad-band SEDs built with simultaneous data taken at 
the epochs of the XMM-Newton observations show a clear UV 
excess. The high UV fluxes are explained if we assume a contri- 
bution by thermal radiation from the accretion disc. 

On the other side, the corresponding X-ray spectra indicate a 
possible soft excess. When comparing our X-ray data with pre- 
vious data from other satellites, the X-ray spectrum appears to 
vary dramatically, so the soft excess cannot be ascribed to the 
accretion disc and/or hot corona surrounding it. A much more 
variable emission contribution is required. 

We found that the broad-band SEDs of BL Lacertae can 
be explained in terms of two synchrotron emission components 
with their corresponding SSC radiation, plus a thermal compo- 
nent representing the contribution of the accretion disc. When 
fitting these two non-thermal components by means of the heli- 
cal jet model of Villata & Raiteril (11999i) . and the thermal one 
with a black body law, we find that the accretion disc has a 
temperature <: 20000 K and a luminosity ^ 6 x lO'^'^ergs 
Taking into account that L = rjMc^, with rj ^ 0.06 in the case of 
Schwarzschild's metric (Shakura & Sunvaev, 1973), we can de- 
rive a lower limit to the accretion rate: M ^ 0.2 Mq yr"'. And if 
we assume that the luminosity equals the Eddington's critical lu- 
minosity, we can also infer a lower limit to the black hole mass: 
Mbh ^ 6 X IO^Mq. This value scales as L^mI^. 

A further emission component coming from inverse- 
Compton scattering on external photons from the disc and/or the 
broad line region might be present, but it is not needed to fit the 
observations, so we think that its possible contribution would be 
a minor one. Its modelling would introduce additional parame- 
ters that we could not reliably constrain. 

Being aware that also other interpretations might be able to 
account for the observations reported in this paper, we notice 
that a helical jet model is motivated by some observing evidence. 
Indeed, VLBA/VLBI studies of the jet structure in AGNs have 
revealed bent jet morphologies that ar e suggestive o f streaming 
motions along a helical path (see e.g. iListeA 12001 1) or that the 
magnetic field may present a helical geometry (Gabuzda et al.', 
2004); this is also true for BL La certae (fTatevama et al., 1998; 
benn etin. l2000VMarscher et al.', 2008; O'Sull ivan & Gabuzdal 
2009 1 see also Sti rling et al.. .2003) . Mo reover, a rotating h eli- 
cal path in a curved jet was invoked by IVillata et al.l (l2009l) to 
explain the optical and radio behaviour of BL Lacertae in the 
last forty years, in particular the alternation of enhanced and 
suppressed optical activity, accompanied by hard and soft radio 
events, respectively. If the jet has a helical structure, different 
portions of the jet may be well aligned with the line of sight, 
with consequent Doppler beaming of the emitted radiation. The 
two emission components proposed in our modelling could cor- 
respond to two of these regions, the higher-energy one being lo- 
cated closer to the jet apex. An analogous picture was suggested 



bv lVillata& Raiteril d 19991) for Mkn 501, whose radio-to-X-ray 
multiepoch SED was explained in terms of two jet regions with 
different curvature. Alternatively, we can imagine that the two 
contributions come from two interweaved helical filaments, but 
we regard this hypothesis as less likely. 

Finally, we mention that also in the case of Mkn 421, a pos- 
sible explanation for its X-ray and optical flux behaviour during 
the June 2008 flare implies the ex istence of two different syn - 
chrotron emitting regions in the jet (iDonnarumma et aTll2009bl) . 

Further multiwavelength observations, including GeV data 
from the Fermi satellite, and TeV data from ground-based 
Cherenkov telescopes, will help verify our interpretation. 
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